Four-Loop Collinear Anomalous Dimension in J\f = 4 Yang-Mills Theory 
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We report a calculation in Af = 4 Yang-Mills of the four-loop term in the collinear anoma- 
lous dimension g(X) which governs the universal subleading infrared structure of gluon scattering 
amplitudes. Using the method of obstructions to extract this quantity from the 1/e singularity in 
the four-gluon iterative relation at four loops, we find g^' = —1240.9 with an estimated numerical 
uncertainty of 0.02%. We also analyze the implication of our result for the strong coupling behavior 
of g(X), finding support for the string theory prediction computed recently by Alday and Maldacena 
using AdS/CFT. 



I. INTRODUCTION 

Gluon scattering amplitudes in QCD and supersym- 
metric gauge theories are notoriously difficult to com- 
pute, yet possess remarkably simple hidden structure. 
We expect the greatest simplicity in the maximally su- 
persymmetric jV = 4 Yang-Mills theory, where planar 
L-loop amplitudes are believed to satisfy iterative rela- 
tions [J, 0, 0, H S, @, 0| in L. These relations allow the 
complete, planar all-loop n-particle MHV amplitude A 
to be written in an exponential form due to Bern, Dixon 
and Smirnov [f|. Their ansatz for the n — 4 particle 
amplitude is 
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where A = gy M N is the 't Hooft coupling, s and t are the 
usual four-particle Mandelstam invariants, and the in- 
frared divergences are encoded in the prefactors Adiv In 
dimensional regularization to D = 4 — 2e the divergences 
take the form 



Adiv{s,t) = cxp 
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in terms of an IR cutoff scale fi (required on dimensional 
grounds) and two functions f^~ 2 \\) and (^"^(A) which 
respectively are related to the /(A) appearing in and 
to a second function called g(X) according to 



dX 



A^ ( - 1} (A). (3) 



The three functions /(A), g(X) and c(A) are universal; 
the same functions appear in the exponential ansatz for 
any planar n-particle MHV amplitude 

The function /(A) is well-known from another role: 
it is the cusp anomalous dimension, which also governs 



the scaling of twist- two operators in the limit of large 

spin s @, la [IE EI, 



A (Tr[ZD s Z}) - S = /(A) log S + O(S 
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At weak coupling it has been computed through four 

loops Em, Hi sis, 
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with r = —2.00002(3) (the quantity in parentheses de- 
notes the current best numerical uncertainty in the last 
digit), while AdS/CFT calculations [IE El indicate the 
strong coupling behavior 
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The recently proposed dressing phase Il8l. lilil ] for the 
asymptotic S-matrix [H H3, HE HE HETll of the spin 
chain description of planar Af = 4 Yang-Mills implies 
that at finite A, /(A) satisfies a certain integral equa- 
tion [26| whose solution is compatible with the lim- 
its ©, © HHH, HI] and moreover predicts that r = -2 
in ©. 

Much less is known about the second function g(X) 
which governs the subleading infrared divergence and 
may be called the "collinear" anomalous dimension [3fj| . 
Perturbative calculations through three loops [f| have es- 
tablished that 
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and the purpose of this note is to report a numerical 
calculation of the four-loop coefficient 
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(4) 



-1240.9(3). 



(8) 



It is an important outstanding problem to relate g(X) to 
more familiar observables which could perhaps be com- 
puted using integrability techniques. In particular it 
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would be interesting to derive an integral equation satis- 
fied by 5(A). 

An important step forward was recently taken by 
Alday and Maldacena [3l|, who gave a prescription 
for computing gluon amplitudes at strong coupling us- 
ing AdS/CFT and found perfect agreement with the 
ansatz ([T]). A consequence of their calculation is the pre- 
diction 



3 (A) = 2(l-log2)^— + 0(1) (9) 

for the leading strong-coupling behavior of g(X). It is 
important to mention that while /(A) is independent of 
the IR renormalization scale, the same is not true for 
5(A). Rescaling jit to kfi sends g(X) to g(X) + 2 log(fc)/(A). 

In section II we review our calculational method, which 
is the same as the one used in [16j to calculate to very 



high accuracy. In section III we explore the connection 
between the weak- and strong-coupling behavior of g(X). 
First we use the Alday-Maldacena result ((9]), together 
with the known data (0 through three loops, to build an 
interpolating function and make a 'prediction' for g^ . 
Remarkably our result © is within 7% of this predicted 
value. Finally we use all currently available data on g(X) 
to construct a Pade approximant which represents our 
best guess for the function's behavior at finite A. 



II. METHOD AND RESULTS 

To compute g^ we apply the method of obstructions 
developed in [l(| [32J. We begin with the four- loop it- 
erative structure encoded in |T]), which in dimensional 
regularization takes the form [5| 



M (4) (e) - M (3) (e) M (1) (e) + M (2) (e) (M (1) (e)f 
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where M^ L \e) is the ratio of the L-loop amplitude to the 
tree amplitude and /q 4 \ fi^ are two numbers. According 
to whose conventions we follow except as noted below, 
these are related to the four-loop coefficients and g^ 
in /(A) and g(X) by 



f( 4 ) = LitW 
Jo 2 4 4 ' 



(ii) 



The factor of 1/2 4 arises here at four loops because the 
expansion parameter used by [j| differs by a factor of 2 
from that of [3l[ which we use in §5§ and |J7J|. 

Each of the five terms on the left-hand side of (fit))) is 
separately an extremely complicated function of the ratio 
x = t/s with singular behavior starting at C(l/e 8 ). It is 
therefore rather remarkable that the five terms conspire 
to add up to such a simple object involving only two free 
constants @, i, Q. The leading singularity of Af (1) (4e) 
on the right-hand side is 0(l/e 2 ), so we can determine 
Z} 4 "* by reading off the coefficient of the l/e pole on both 
sides of (|TU|. 

The idea behind the method of obstructions, explained 
in detail in [16), is that we do not need to fully compute 
all of the M^ L \e). Rather, it is sufficient to compute 
what are in some sense the 'constant pieces', since the x- 
dependent pieces are guaranteed to cancel each other out 
on the left-hand side of (fTI)]) . This idea is made precise 
by writing each amplitude as a Mellin transform, in the 
form 
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for some F(y, e). It turns out that f[ 4 ^ multiplies 6(y) in 
Mellin space, so in order to isolate this term it is sufficient 
to read off just the coefficient of S(y) on both sides of (|10p . 
throwing away the rest. We call terms proportional to 
(5-functions obstructions. The utility of this method ben- 
efits significantly from the fact that we can work directly 
in Mellin space, where it is relatively easy to construct 
explicit formulas for Feynman integrals [H, [34| • 

It was further shown in [l|| that obstructions obey a 
product algebra structure, meaning that the obstruction 
in any product of amplitudes is given by the product of 
their obstructions. This statement is obvious in y-space 
where it hinges on the simple fact that S(y) convolved 
with itself is again 5(y). Consequently we can calculate 
the obstruction p( L \e) in each M^(e) separately and 
then insert those results into the polynomial in (JTUJ) . 

The L — 1,2,3,4 loop amplitudes we require may be 
expressed in terms of the 12 scalar Feynman diagrams 
depicted in Fig. 1 I5I Il5l. l35l . l36j . The reader may find all 
necessary details in [g, ll5( . Our convention is that each 
loop momentum integral comes with a factor 
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The standard convention includes only the factor in 
brackets. The additional factor (— s) e / 2 (— t) e / 2 renders 
all amplitudes dimensionless (when the appropriate nu- 
merator factors are included), but it does not alter the 
form of (|10[) since each term in that equation has a com- 
mon factor (— s)~ 2e (— t)~ 2e . The factor of —1/2 in front 
here eliminates the need for a factor of (— 1/2) L in H 
below. 
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FIG. 1: The 12 integrals appearing in the L < 4-loop four- 
particle amplitudes in J\f = 4 Yang-Mills theory. We refer the 
reader to [15|] for all necessary details. 



We have computed the obstructions in the 12 separate 
integrals through the first eight orders in e and present 
the results in equation (fTT)]) at the end of the paper. The 
calculation was performed using the algorithm explained 
in [l6[ based on Czakon's MB program [37|. Numerical 
integrations were performed using CUBA's Cuhre algo- 
rithm ,38], and the digit in parentheses denotes the re- 
ported uncertainty in the final digit. 

An exact formula for 2"W(e) was given in [l6j], as were 
the first seven terms in the two- and three-loop ampli- 
tudes (though here we indicate separately the contribu- 
tions from the two three- loop integrals). No significant 
effort was spent on attempting to improve the numerical 
accuracy of the C(l/e 3 ) and 0(l/e 2 ) terms in the indi- 
vidual four-loop integrals because these values play no 
role in this paper. 

To find the full obstruction p( L ) in the L-loop ampli- 
tude we add together the contributions from the individ- 
ual integrals according to d, [H, HE Hll 



p(i) = i« 

p(2) = 2I (2) ; 

p( 3 ) — 21^ a + 4j( 3)b 

p(4) = 2T (A)a +42 -(4)6' +4 j(4)c 



+■ 21^° 

.j(4)/2_ 



8I (4)e 



(14) 



Plugging the resulting expressions for P( L ) into (HO]) m 
place of M^ L ' and using the relation (fTTj) leads to the 
advertised value ([5]). 



III. BRIDGING WEAK AND STRONG 
COUPLING 

As mentioned above, string theory provides predictions 
for how /(A) and g{X) behave at strong coupling. It is 
therefore tempting to try to make a connection between 
the weak and strong coupling regimes. An approxima- 
tion scheme that has been very successful for /(A) is the 
following. If n + 1 orders of the expansion around A = 
are known, then one constructs an approximating func- 
tion /(A) as an appropriate solution to the polynomial 
equation (3f| 
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16tt 2 
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c r f(\y, 
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where the coefficients c r are determined by imposing that 
/(A) agree with /(A) through 0(\ n+1 ). 

The form of (fl~5)) incorporates a strong coupling expan- 
sion for /(A) in the parameter 47r/\/A with leading term 
aVX/4ir where a 2n — l/c 2n . 

Using the expansion of /(A) up to four loops one finds 
an approximating function /(A) that agrees, within a few 
percent, with the string theory strong coupling prediction 
and with the BES ansatz |26j for all positive values of the 
coupling. 

The success of this approximation scheme makes its 
application to g(X) a natural thing to attempt. How- 
ever, since g(X) is negative for small A but becomes pos- 
itive at strong coupling, it must have at least one zero 
for A > 0. Clearly, any approximating function of the 
form (|15|) (with an appropriate modification because 5(A) 
only starts at two loops) 
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(16) 



leads to a g(X) that cannot possibly have a zero for A 7^ 
and hence a contradiction. One way to find a consistent 
approximating function was introduced in 31]. The idea 
is to perform a change of the IR renormalization scale 
li — > exp(£/2)/i. Then, as discussed above, g(X) — * g(A) + 
£/(A). The new function g(X, £) = g(X) +£/(A) can now 
be approximated by using (|15p for any £ > 0. 

Consider the approximation (|15p for <?(A, £) with n = 2, 
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16tt 



= c 2 f (A, 2 + c 3 |(A, 3 + c 4 ^(A, 4 . (17) 



This equation depends on four parameters, namely, 
C2)C3,C4 and £. Imposing that the approximating func- 
tion g(X, £) agrees with g(X, £) through C(A 3 ) determines 
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able u = \/\ + X/tt 2 was considered for /(A). Specifically 
we consider the ansatz 



FIG. 2: A Pade approximant for the function g(X) versus 
A/I6-7T 2 , obtained by fitting the ansatz (|I8|) to all available 
data at weak and strong coupling. 



the three coefficients 02,03,04 in terms of £. In [31[ 
the resulting approximating function was extrapolated to 
strong coupling and compared to the string theory pre- 
diction for a range of £. At the special value £ = \ log 2, 
for example, the approximation gives g(X) — 1.37\/A/47r 
compared to the string theory prediction g(X) — 2^/X/Att. 

Here we would like to use the results of Alday and 
Maldacena to make a prediction for g^ to compare our 
result ([8]) to. Imposing that g agrees with the string 
theory prediction at strong coupling fixes the parameter 
£ « 0.73679. Then all four parameters 02, C3, C4 and £ are 
completely fixed and we can expand the resulting g in A 
to C(A 4 ) to find the predicted value « -1336.9. Re- 
markably our result 1(5)1 is only 7% away from this value. 
We interpret this as good evidence for the string th eory 
strong coupling prediction of Alday and Maldacena [3lj . 

Finally we present also a more conventional Pade ap- 
proximant for .g(A) based upon all currently available 
data. Here we follow the approach described in [l5[ where 
a [3/2] Pade approximant in terms of the auxiliary vari- 



G(A) 



N + Niu 
1 + D x u + D 2 u 2 



(18) 



with four parameters iVo, iVi, D%, D%- The choice u = 
was motivated by evidence that /(A) should 
have a branch point at A = —it 2 . It seems reasonable to 
guess that the same might be true for 5(A), though we 
have no direct evidence for this guess. 

The four parameters in (|18p can be uniquely deter- 
mined by fitting to all available data — the perturba- 
tive expansion through four loops as well as the strong- 
coupling limit ([9]). The resulting approximant, displayed 
in Fig. 2, may be considered our best candidate picture 
of 5(A) at the moment. Note that it is minimal in the 
sense that G(A) was designed to have a single zero along 
the positive real axis, whereas the true 17(A) could poten- 
tially have any odd number of zeros. Curiously, the zero 
lies very close to A/167T 2 = 1, although it is impossible 
to conclude based on the limited available data whether 
or not this is just a coincidence. 
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